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Dynamics of unvisited sites in presence of mutually repulsive random walkers 
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We have considered the persistence of unvisited sites of a lattice, i.e., the probabihty S{t) that a 
site remains unvisited till time t in presence of mutually repulsive random walkers. The dynamics of 
this system has direct correspondence to that of the domain walls in a certain system of Ising spins 
where the number of domain walls become fixed following a zero termperature quench. Here we 
get the result that S{t) oc exp{~at^) where f3 is close to 0.5 and a a function of the density of the 
walkers p. The number of persistent sites in presence of independent walkers of density p' is known 

to be S'(t) = exp{—2'^J^p't^^^). We show that a mapping of the interacting walkers' problem to 

the independent walkers' problem is possible with p' = p/(l — p) provided p',p are small. We also 
discuss some other intricate results obtained in the interacting walkers' case. 

PACS numbers: 05.40.Fb,05.50.+q,02.50-r 



I. INTRODUCTION: THE ORIGINAL SPIN 
PROBLEM 

Dynamical evolution of a spin system following a 
quench to zero temperature from a disordered state may 
lead to a non-equilibrium state, e.g., as in the one di- 
mensional ANNNI (Axial Next Nearest Neighbor Ising) 
model [H which has the Hamiltonian 



(1) 



where 5*^ = ±1 is the spin at the ith site and k > 
is the ratio of the second neighbour and first neighbour 
interactions. For k < 1, the quench does not lead to 
the equilibrium configuration which is ferromagnetic for 
K < 0.5 and antiphase for k > 0.5. Starting from a 
random state, there is a short initial time during which 
the domains of size one die and this eventually results in 
a configuration with fixed number of domain walls. In 
this state, the domain walls become "fluid" in the sense 
that they can move indefinitely (keeping the energy of 
the system same) 0, 13, 01 but cannot cross each other. 
In such a system the persistence dynamics shows that 
the number of persistent spins is neither a power law nor 
exponential but rather follows a stretched exponential 
decay, 
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(2) 



with a « 1 and /3 = 0.45 |3i]. 

The above dynamical scenario can easily be repre- 
sented by an equivalent system of mutually avoiding ran- 
dom walkers and the fraction of persistent spins will then 
be given by the fraction of unvisited sites SYiVtill time t 
in the system g i 0, S i, M MMMM,MM^- 
Representation of the spin dynamics by random walk of 
domain walls is well-known; in the common examples like 
the Ising or Potts model, the random walk is accompa- 
nied by annihilation of two domain walls if they meet 



In the original ANNNI model problem, when there are 
M domains with each separated by at least two lattice 
spacings, the probability distribution of the number of 
domain walls M within a size L is 



P{L,M) 
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This equation is easy to derive once it is realised that the 
problem is identical to the Bose statistics of distributing 
L particles in M boxes with the number of particles in 
each box greater than or equal to 2. It has been observed 
that the ratio po = M/L has a mean value quite close 
to the most probable value po — 0.2764 Distribution 
of the value of po can be calculated numerically which 
shows that its fluctuation decreases with the system size 
(typically Apo decreases from 0.0645 for L = 20 to 0.0225 
for L — 175 in a power law manner, Apo ^ iV^*^'^). This 
indicates that the ANNNI dynamics reflects the behav- 
ior of S{t) for a specific value of the density of walkers 
in the equivalent random wall picture. It is therefore a 
meaningful exercise to find out S(t) for general density 
p = N/L where N is the number of mutually exclusive 
walkers and L the chain length. In this paper we have 
considered < p < 1 and calculated numerically S{t). 
S(t) shows a stretched exponential behaviour with an 
exponent close to 1/2. 

This problem, as we have shown, can be mapped to 
that of the independent random walkers with a subtle 
difference and also for small p. The latter problem has 
been exactly solved [isj and we compare its results with 
that of the numerical simulation of mutually repulsive 
walkers to find a good agreement at small values of p. 

Our results also show a difference when persistence is 
calculated in terms of a spin system and the system of 
pure brownian walkers which we have discussed in the 
paper. Some other dynamical properties have been in- 
vestigated in this context. 
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FIG. 1: Possible movement of the tagged walker (highlighted): 
(a) cannot move, (b) either does not move or can move to the 
left/right, (c) either does not move or can move to the right, 
(d) either does not move or can move to the left. 



II. PERSISTENCE OF UNVISITED SITES IN 
PRESENCE OF MUTUALLY EXCLUSIVE 
WALKERS 

To keep correspondence with the original spin dynam- 
ics of the ANNNI model, any of the site in the system 
may be selected for updating. It may or may not contain 
a random walker. If it does, then the step taken by this 
random walker is according to the following three situa- 
tions: 

1 . If the walker is flanked by two random walkers on both 
sides; it cannot make any movement and stays there - this 
corresponds to a locked domain wall (Fig la). 

2. There are no neighbouring random walker, then the 
random walker remains at its position with probability 
1/2 - this corresponds to the probability that a spin does 
not flip even at the domain boundary. It can also move 
to either left or right with equal probability (Fig lb). 

3. There is only one neighbouring walker, say, to the left 
(right), then it does not move or moves to the right (left) 
(Figl (c) and (d)) with equal probability. 

Let us represent this dynamical rule by Dann, a dynam- 
ics which corresponds to the ANNNI model dynamics. 

As the walls are reflecting, if the random walker hap- 
pens to hit the wall, it can only step in-wards or stay 
there. When L sites are hit, one Monte Carlo step is said 
to be completed. 



It may be noted that in the original ANNNI model, 
domain walls need to maintain a least distance of two 
lattice spacings. We have relaxed this condition to one 
lattice spacing which is equivalent to having on-site re- 
pulsion of random walkers. Effectively, this means that 
the original ratio po — M/L in ANNNI corresponds to 
p — N / L = 2/3o in the present case. 

The starting position of the random walkers may be 
assumed to be either already visited or not yet visited. 
The calculation of persistence will depend on it. In what 
we call the spin picture (SP), they are unvisited and in 
the random walker picture (RWP) they are assumed to 



be visited already. This brings in a subtle difference in 
the two problems. We have discussed the two problems 
separately in the following two subsections. 

A. The spin picture 

First, we calculate the survival probability S{t) defined 
as the probability that a site has not been visited by any 
of the random walkers till time t. (This corresponds to 
P(t), the persistence probability of the original ANNNI 
model with p ~ 0.54.) 

In Fig. 2, S{t) is plotted against time t for different 
densities p with a fixed lattice size L — 10000. In each 
case, S{t) decays with time following a stretched expo- 
nential behavior exp(— af^) where (i — 0.50 ± 0.02. This 
value of P compares well with 0.45 obtained for P{t) in 
the ANNNI model [3]. Increasing the system size does 
not affect the result, only the variation of S{t) can be 
observed over a longer duration of time. 



Interestingly, the value of a increases with the the in- 
crease of the density p till the value of p w 0.55. Beyond 
this value, a decreases gradually as the density increases 
roughly as (1 — p) (Fig. 3). The behaviour of a as a func- 
tion of p will be discussed in greater detail in the next 
section. 

The qualitative behaviour of a as a function of p is not 
difficult to explain. For small values of p, the probability 
that a domain wall is hit is small and therefore one gets 
a slow decay of S{t) with time. On the other hand, for 
large values of p, most of the domain walls will be 'locked' 
such that again the variation of S{t) will be slow. 

For very small p, we have checked that 

ail -p) = V2a(p) (4) 
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FIG. 2: (a) Survival probability S{t) as a function of 
time t for different densities of walkers on a lattice of size 
10000. Each curve follows a stretched exponential of the form 
exp(-at°-^''). 
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FIG. 3: Slope of the stretched exponential curve a is plotted 
against the density p for both the SP (spin picture) and RWP 
(random walkers picture) for the dynamical rule Dann (see 
section II). Both the curves follow a behaviour 1.13p/(l — p) 
for small p, shown by the dashed line. 

holds good to a high degree of accuracy. Effectively this 
means that the relaxation rate at p ^ 1 is twice as much 
as that at p — > 0. This can be explained in the following 
manner: At p — > 1, an empty site occurs with a probabil- 
ity (1 — p) but has both neighbours occupied with a very 
high probability, almost one. At p ^ 0, an empty site 
having an occupied neighbour is very rare (probability 
proportional to p). In the former, in one iteration, the 
empty site can be visited by either of its neighbouring 
walkers while in the latter, visit to the empty site is pos- 
sible by one walker only, making the time scale double. 



One can now compare the value of a obtained in the 
present study with that of the ANNNI model where p ~ 
0.27. In [sl, the value of a was found to be equal to 1.06. 
For p = 2po ~ 0.54, we find that a is very close to the 
value 1 (Fig. 3) showing again a good agreement with 
the ANNNI model dynamics. 

B. The Random walker picture 

In the RWP everything remains same but the ini- 
tial sites occupied by the random walkers are as- 
sumed to be already visited. Now we find that S{t) 
again has a stretched exponential behaviour : S'(i) ~ 
exp{—aRwt'^'^^) with Prw = 0.5±0.01. Now ajiw does 
not show any non-monotonic behaviour but appears to 
diverge as p ^ 1. This is again understandable, in the 
present picture, when p is close to one, most of the sites 
are already non-persistent to begin with and S(t) decays 
very fast making urw oo- 

The exponents (Srw and /3 are apparently equal in the 
two pictures. In Fig. 3 the behaviour of aaw against p 
is also shown. It is to be noted that upto p w 0.5, a and 



c^RW are equal and behave differently beyond this point. 
For the SP, the possibility of the domains getting locked 
increase as p increases and this happens with a higher 
probability beyond p = 0.5. 

III. MAPPING TO A SYSTEM OF 
NON-INTERACTING WALKERS 

In the last section we obtained the result that the frac- 
tion of unvisited sites S(t) in presence of mutually repul- 
sive walkers has a stretched exponential decay with expo- 
nent 1/2 in both the SP and RWP. This behaviour turns 
out to be exactly the same as that of S'{t), the num- 
ber of unvisited sites in the presence of independent or 
non-interacting walkers. In the latter system it has been 
shown pJi] that when p' is the density of independent 
walkers, 

S'{t) - exp{~a'p'Vt), (5) 



with a' = 2 




In this section we show that the interacting system 
can be mapped to the independent walkers' system with 
the transformation p' = p/{l — p). To show this, let us 
consider a configuration C of interacting walkers of 
density p which follow a dynamics represented by D. For 
the present discussion, we make the dynamics D simpler 
than Dann- the walker will always execute a movement 
when at least one of the neighbouring site is vacant - if 
both are vacant, probability to move either to the left or 
to the right is 1/2. In case only one neighbouring site 
is occupied, the random walker will move to the empty 
neighbouring site. (One obtains the same behaviour of 
S{t) with this rule, including the relation only the 
numerical value of a increases by a factor of \/2 compared 
to Dan n where the time scale is simply double compared 
to D.) 

For the independent walkers, let us consider a config- 
uration Co of walkers of density p' , who do not "see" 
each other. The dynamics Dq here is simply that each 
walker will move to left or right with equal probability. 
The world lines in the 1-1-1 dimension of the walkers are 
shown in Fig 4a and 4b. 

Now let us create a mapping of the original configura- 
tion C to C" given by 

x'^{t) ^ Xkit) - k, (6) 

where Xk{t) is the position of the fcth walker (fc — 
1,2,...M from the left) at time t Effectively this 

mapping implies that one spacing between consecutive 
walkers is being removed. This would remove the con- 
straint in C that two walkers have hard core repul- 
sion and each world line of C therefore also occurs in 
Co. Although all world lines of C and Co have one to 
one correspondence, in C' one has the constraint that 
xi < X2 < X3 < xn while in Cq there is no such con- 
straint. Therefore a particular configuration may occur 
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FIG. 4: Typical examples of (a) mutually repulsive walkers 
(C) and (b) Independent walkers (Co), (c) The mapped con- 
figuration C" following eq. 6. In (d), another realisation for 
the independent walkers is possible by switching the order of 
the walkers at the later time. This is not allowed in C". 

with different weight factors in Co and C (Figures 4b, c, 
and d). 



Obviously a "three" state cannot be avoided if p > 2/3 

and this gives an upper bound where the disagreement 
will occur. In reality, such states occur at values of p 
much below than this, even at about p = 0.2. We have 
verified that, if the three-states are forcibly ruled out in 
the simulation, the correspondence between the indepen- 
dent and interacting walkers remain valid upto p « 0.4. 

We would like to comment in this section that while 
for the interacting walkers' system, a behaves differently 
in the SP and RWP, no such difference exists for the 
independent walkers' case. This is because there is no 
restriction on the movement of the walkers here even as 
the density becomes high. 

A subtle point relevant to the mapping needs to be 
mentioned here. At small p the results for C and Co 
are equivalent indicating that tlic^ dynamical evolution of 
the walkers can be mapped to each other. It may still 
remain a question whether the persistence probability of 
C can be mapped to that of C. The question arises as 
the N sites removed from the original system may either 
be persistent or non-persistent. On an average, however, 
the persistence of the two systems C and C will be same. 
This is because the average number of persistent sites re- 
moved is P{t)N. Thus in the mapped system, persistence 
probabilty is again (P(t)L - P{t)N)/{L - N) ^ P{t). 
This justifies the correspondence of persistence in C and 
C and hence Cq. The issue of equivalence of persistence 
requires special mention as persistence is not related to 
other dynamical behaviour of a system in general. 



Using the above mapping, the effective chain length in 
C is L — N and not L. Thus the density p in C is related 
to the density p/ in C by the equation 



Ignoring the weight factor, the mapping is effective in 
showing the correspondence between the interacting and 
independent walkers' picture. Exact correspondence will 
imply that a or auw (with D) would be equal to 

where a' = 2y^(2/w), when dynamics D is used. This can 
happen if the dynamical rule D applied to C leads to 
states which when mapped to C will correspond exactly 
to the states obtained by applying Dq on Co (with the 
same weightage). We have verified that this is true for 
configurations in which either a walker is "alone" (both 
neighbours are empty) or has at most one walker in a 
neighbouring site. However, when three walkers occupy 
consecutive sites (a "three" state), the dynamics D gives 
rise to states which cannot be obtained from Cq applying 
Do on it. Since the probability of having "three" states 
increases with p, wc expect that the results for indepen- 
dent and interacting walkers will differ at higher p. For 
Dann, it is expected that a and a^w values would be 
equal to ■^p/(l — p) = upto small p which is ex- 

actly what we observe (Fig. 3) (time scales for Dann 
being simply twice that of D). 



IV. NON-MONOTONICITY OF a AND A FEW 
RELEVANT COMMENTS 

The result for persistence probability in the spin pic- 
ture and random walker picture differ in the interacting 
walkers' case as in the SP there is a non-monotonicity in 
a. This non-monotonic behaviour is clearly due to two 
features (a) presence of interacting walkers and (b) the 
dynamic quantity under consideration being persistence. 

Point (a) is already discussed in the last section. 
Regarding point (b), it must be noted that the non- 
monotonicity appears when we assume that the initially 
occupied points are not visited, a fact which is relevant 
to persistence dynamics only. In this section we have dis- 
cussed a few other dynamical phenomena in presence of 
interacting walkers. However, we find that none of these 
are accompanied by any non-monotonic behaviour of the 
relevant quantities appearing in them. 

Two dynamic quantities ai and (J2 representing fluc- 
tuations can be defined in the following way: we tag a 
random walker and calculate the fluctuation of its posi- 
tion x{t) at time t with respect to its initial position x{0) 
and study its behavior with time (Fig. 5). 

Precisely, a\ is defined as 

cj,{t) = ^{{<t)-x{Q)f). (8) 
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FIG. 5: The movement of a tagged domain wall along time 
(vertical axis): a;(0) is its initial position, x{t) its position at 
time t and {x{t)), the mean position averaged up to time t. 



In the second measure, we notice that the path of a walker 
can be viewed as an interface (with no overhangs) in 1 + 
1 dimensions (Fig 5). One can then measure the interface 
width at any time given by 



^2 



(9) 



where {x{t)) is the mean value of the position x at time 
t. It is known that for a single walker (i.e., p = in the 
thermodynamic limit) ai{t) = At^ with 9 ~ 0.5. Here, 
we find o'i(t) — with 9 ~ 0.25 at longtimes for all 
values of p. This is in agreement with [20| where the 
result 9 = 0.25 has been derived exactly. In the present 
system, p has a finite value and for the smallest value of 
p shown in Fig. 6, a crossover effect is noted, i.e., the 
behavior at earlier time appears to be consistent with 
t"'^. This is because at small p, the walker continues as a 
free walker for a considerable period of time and exhibits 
the corresponding behaviour. 



The behavior of (Ti{t) with time t has been studied for 
values of p even smaller than 0.1 in smaller lattices and 
it appears that for any finite p, however small, 9 ~ 0.25 
is valid at longer times always. We conclude that there 
is a transition point at p = for any p ^ 0, the random 
walker exponent is ~ 0.25 while for p = 0, it is 0.5. The 
results for <72 are consistent with the above observations. 
Both CTi, (72 — At^ with 9 = 0.25 independent of p (for 
p ^ 0), while A depends on p. In Fig. 7, we plot A{p) 
against p for both cri and a2- A{p) decereases monotoni- 
cally with p and follows a rough exponential decrease as 
A{p) ~ exp(— 2p) except for values very close to 1, where 
one can expect anomalous behavior. 

We investigate the behavior of another quantity D{t), 
which we define as the fluctuation of the distance d{t) 
between two neighbouring walkers at time t with respect 
to its initial value d{0). Precisely, 
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FIG. 6: Fluctuation ai of the position of a walker w.r.t. its 
initial position as a function of time t for different no. of 
walkers on a lattice of size 5000. The best fit curves have a 
slope ~ 0.25 for the higher densities. 
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P 

FIG. 7: Amplitudes (A) of ai,a2{— At^) and Dsat as function 
of p. The dashed lines are guides to the eye. 



D{t) shows an initial increase with t and reaches a time 
independent equilibrium value D^at at larger times. This 
equilibrium value D^at (calculated from the mean value 
of the last 100 Monte Carlo steps), when plotted against 
p, again shows a monotonic decay with p (Fig. 7). Hence, 
in contrast to the factor a appearing in the persistence 
dynamics, we do not find non- monotonic behavior in the 
factors appearing in the other dynamical features. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have considered the dynamics of N 
mutually avoiding random walkers on a one dimensional 
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chain of length L which simulates the quenching dynam- 
ics in the ANNNI model for a particular value of p ~ 0.54. 
For this value of p, we verify that the survival probability 
S{t) which corresponds to the persistence probability in 
the ANNNI model follows a stretched exponential behav- 
ior consistent with the ANNNI model dynamics. Wc also 
observe very good quantitative agreement for the expo- 
nent (3 and slope a. On generalising the value of p, we 
find that the behavior S{t) ~ exp(— ai'') is valid for all 
p with (3 showing a universal value of 0.50 ± .02. 

Observing that the time dependence in S{t) is identical 
to that in S'{t) (the corresponding quantity in presence 
of non-interacting walkers) for all p, we have shown that 
a mapping between the two indeed exist which remains 
exact for small p as far as the behaviour of a is concerned. 

We have considered two different pictures while com- 
puting the persistence probability; in the spin (random 
walker) picture the sites initially occupied by the walkers 
are assumed to be not visited (visited) and the behaviour 
of a as a function of p is sensitive to this difference. In 
the SP, it has a non-monotonic behaviour. Such non- 
monotonic behaviour emerges in the case of interacting 



walkers only. However, when other dynamical phenom- 
ena in presence of interacting walkers are studied, no such 
non-monotonic behaviour is found. Thus we find that the 
persistence dynamics in a system with a finite density of 
mutually avoiding random walkers calculated in terms of 
the fraction of sites unvisited till time i, shows a unique 
behavior compared to other dynamical quantities. This 
again supports the fact that persistence is a phenomenon 
which cannot be directly connected to other dynamical 
features of a system. 
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